Kinematic dynamos using constrained transport with high order Godunov schemes and adaptive mesh refinement
نویسندگان
چکیده
We propose to extend the well-known MUSCL-Hancock scheme for Euler equations to the induction equation modeling the magnetic field evolution in kinematic dynamo problems. The scheme is based on an integral form of the underlying conservation law which, in our formulation, results in a “finite-surface” scheme for the induction equation. This naturally leads to the well-known “constrained transport” method, with additional continuity requirement on the magnetic field representation. The second ingredient in the MUSCL scheme is the predictor step that ensures second order accuracy both in space and time. We explore specific constraints that the mathematical properties of the induction equations place on this predictor step, showing that three possible variants can be considered. We show that the most aggressive formulations (referred to as C-MUSCL and U-MUSCL) reach the same level of accuracy as the other one (referred to as Runge-Kutta), at a lower computational cost. More interestingly, these two schemes are compatible with the Adaptive Mesh Refinement (AMR) framework. It has been implemented in the AMR code RAMSES. It offers a novel and efficient implementation of a second order scheme for the induction equation. We have tested it by solving two kinematic dynamo problems in the low diffusion limit. The construction of this scheme for the induction equation constitutes a step towards solving the full MHD set of equations using an extension of our current methodology.
منابع مشابه
A High Order Godunov Scheme with Constrained Transport and Adaptive Mesh Refinement for Astrophysical MHD
Aims. In this paper, we present a new method to perform numerical simulations of astrophysical MHD flows using the Adaptive Mesh Refinement framework and Constrained Transport. Methods. The algorithm is based on a previous work in which the MUSCL–Hancock scheme was used to evolve the induction equation. In this paper, we detail the extension of this scheme to the full MHD equations and discuss ...
متن کاملA Stable, Accurate Methodology for High Mach Number, Strong Magnetic Field Mhd Turbulence with Adaptive Mesh Refinement: Resolution and Refinement Studies
Performing a stable, long duration simulation of driven MHD turbulence with a high thermal Mach number and a strong initial magnetic field is a challenge to high-order Godunov ideal MHD schemes because of the difficulty in guaranteeing positivity of the density and pressure. We have implemented a robust combination of reconstruction schemes, Riemann solvers, limiters, and Constrained Transport ...
متن کاملNew adaptive interpolation schemes for efficient meshbased motion estimation
Motion estimation and compensation is an essential part of existing video coding systems. The mesh-based motion estimation (MME) produces smoother motion field, better subjective quality (free from blocking artifacts), and higher peak signal-to-noise ratio (PSNR) in many cases, especially at low bitrate video communications, compared to the conventional block matching algorithm (BMA). Howev...
متن کاملAn Adaptive, Higher-order Godunov Method for Gas Dynamics in Three-dimensional Orthogonal Curvilinear Coordinates
We describe an adaptive higher-order Godunov method for the compressible Euler equations with an arbitrary convex equation of state in a three-dimensional system of orthogonal curvilin-ear coordinates. The single grid algorithm is a fractional-step method which uses a second-order Godunov method for gas dynamics in each fractional step. The single grid algorithm is coupled to a conservative loc...
متن کاملar X iv : a st ro - p h / 98 07 12 1 v 1 1 3 Ju l 1 99 8 COSMOLOGICAL ADAPTIVE MESH REFINEMENT
We describe a grid-based numerical method for 3D hydrody-namic cosmological simulations which is adaptive in space and time and combines the best features of higher order–accurate Godunov schemes for Eulerian hydrodynamics with adaptive particle–mesh methods for collision-less particles. The basis for our method is the structured adaptive mesh refinement (AMR) algorithm of Berger & Collela (198...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید
ثبت ناماگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید
ورودعنوان ژورنال:
- J. Comput. Physics
دوره 218 شماره
صفحات -
تاریخ انتشار 2006